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SPECIFICATION 

DATA INTERPOLATING DEVICE AND METHOD, SAMPLING FUNCTION 
GENERATING DEVICE, DATA INTERPOLATING PROGRAM, AND RECORDED 
MEDIUM 

[Background of the Invention] 
[Field of the Invention] 

The present invention relates to a data interpolating 
device and method, a sampling function generating device, a 
data interpolating program and a recorded medium, and more 
particularly to a data interpolating device and method, a 
sampling function generating device, a data interpolating 
program and a recorded medium that are suitable for 
interpolation of the discrete data such as voice or image. 

[Description of the Related Art] 

Conventionally, various data interpolating methods for 
interpolating a value between given sampling values have been 
proposed. The simplest method is a linear interpolation. In 
this linear interpolation, interpolation data represented by 
(xl+x2)/2 is employed as an intermediate value interpolating 
between two sampled values xl and x2, for example. 

Also, a data interpolating method employing a given 
sampling function is well known. For the data interpolation 
in this case, the sampling function as referred to as a sine 
function is conventionally employed. 
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FIG. 1 is an explanatory view of a sine function. The 
sine function appears when a Dirac's delta function is 
subjected to inverse Fourier transform, and is defined such 
as sin(jtf t)/(jrft) where the sampling frequency is f . This 
sine function has the value 1 only at the sampling point t=0, 
and value 0 at all other sampling points. 

FIG. 2 is an explanatory view of data interpolation using 
a sampling function as shown in FIG. 1. In FIG. 2, consider 
an instance wherein the interpolated value y (tO) corresponding 
to a predetermined position tO (distance a from t2) between 
sampling points t2 and t3, when the values of discrete data 
at sampling points tl, t2, t3 and t4 of equal interval are 
Y(tl), Y(t2), Y(t3) and Y(t4), respectively. 

Generally, to determine the interpolated value y(tO) 
using the sampling function, the values (indicated by x sign) 
at the interpolation position to for a plurality of sampling 
functions having the amplitude corresponding to the discrete 
data values Y(tl) to Y(t4) at the sampling points tl to t4 
are determined, and added by convolution operation. 

This data interpolation is employed for a digital-analog 
conversion process for converting discrete digital data to 
a continuous analog signal, and an expansion process for 
expanding the data discretized at the time of compression to 
continuous original data by interpolation. In this case, it 
is required to interpolate the value between given discrete 
sampling points as smoothly as possible. 
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To deal with this requirement, in a conventional data 
interpolating method for artificially increasing the sampling 
frequency through the interpolation process , a stepwise signal 
waveform was generated by holding the interpolated values in 
a sample and hold circuit , and passed through a low pass filter 
to output a smooth signal. 

However, this method had a problem that the output 
continuous signal is degraded in phase characteristic owing 
to the low pass filter. 

Since the sine function converges to 0 as t goes to ±oc , 
to acquire the correct interpolated values, it is required 
to add the values of the sine function corresponding to all 
the discrete data values. In practice, the convolution 
operation is made in a limited range of discrete data considered 
owing to the processing capability and circuit scale. 
Therefore , there was a problem that the determined interpolated 
value contained a truncation error, and the correct 
interpolated values were not determined. 

To solve these problems, the present applicant has 
invented a data interpolating system for providing the output 
waveform with less distortion without degraded phase 
characteristic or truncation error caused by the low pass 
filter, and already filed a patent ( Japanese Patent Application 
No. 11-173245) . The invention as described in Japanese Patent 
Application No. 11-173245 involves processing a digital basic 
function of FIG. 15 to have the amplitudes corresponding to 
n pieces of input discrete data, and synthesizing them with 
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the over-sampling and convolution operation, and as a result, 
making data interpolation based on the sampling function as 
shown in FIG. 3. 

The sampling function as shown in FIG. 3 is a function 
of finite supports dif f erentiable once over the whole range, 
in which the sampling position t along the transverse axis 
has a finite value other than zero only within a finite local 
area, and all zero in the other area. 

Also, the sampling function of FIG. 3 has a feature that 
the maximal value is at sampling point 13 alone, and the value 
is zero at four sampling points tl, t2, t4 and t5, and passes 
through all the sampling points required to acquire the smooth 
waveform signal. 

In this manner , the function as shown in FIG . 3 is a sampling 
function of finite supports that is dif f erentiable once over 
the whole range, and converges to zero at sampling points tl 
and t5. Accordingly, employing the sampling function of FIG. 
3, instead of the conventional sine function as shown in FIG. 
1, to make superposition of each discrete data, the value 
between discrete data is smoothly interpolated using a first 
order dif f erentiable function. Thereby, the low pass filter 
is dispensed with, and the degradation in phase characteristic 
is prevented. 

Since the sampling function of FIG. 3 converges to zero 
at sampling points tl and t5, it is only necessary to take 
into consideration only the discrete data in a range from tl 
to t5 . Accordingly, when one interpolated value is determined. 
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a limited number of discrete data are only considered, thereby 
greatly reducing the processing amount. And though the 
discrete data outside the range from 1 1 to 1 5 should be naturally 
considered, but are not ignored in consideration of the 
processing amount or precision. Theoretically, there is no 
need for considering them, whereby it is possible to prevent 
a truncation error from occurring. 

However, the data interpolating system as described in 
Japanese Patent Application No. 11-173245 provides the 
excellent characteristics as described above , but was required 
to make more complex calculations with the convolution 
operation or moving average operation than the simple data 
interpolating method such as linear interpolation. 

Thus, it is an object of the present invention to provide 
a smooth data interpolation based on a sampling function of 
finite supports dif f erentiable once or more times over the 
whole range through the simpler calculation, thereby 
shortening the calculation time and simplifying the circuit 
for making the calculation. 

[Summary of the Invention] 

The present invention provides a data interpolating 
device comprising plural stages of delay circuits for 
successively delaying discrete data sequentially inputted, 
and a multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the plural 
stages of delay circuits according to the value of a digital 
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basic function and thereby determines interpolation data for 
the discrete data sequentially inputted. 

In the data interpolating device, a plurality of 
over- sampling circuits may be cascade connected, when one set 
of over- sampling circuit consists of the plural stages of delay 
circuits and the multiplication/addition circuit. 

Also , this invention provides a data interpolating device 
comprising three stages of delay circuits for successively 
delaying discrete data sequentially inputted, and a 
multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the three 
stages of delay circuits according to the value of a digital 
basic function and thereby determines interpolation data for 
the discrete data sequentially inputted. 

In the data interpolating device, a plurality of 
over- sampling circuits may be cascade connected, when one set 
of over-sampling circuit consists of the three stages of delay 
circuits and the multiplication/addition circuit. 

Also, the invention provides a data interpolating device 
comprising plural stages of delay circuits for successively 
delaying discrete data sequentially inputted, a 
multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the plural 
stages of delay circuits according to the value of a digital 
basic function and thereby determines interpolation data for 
the discrete data sequentially inputted, and an averaging 
circuit for calculating the average data of adjacent 



interpolation data output from the multiplication/addition 
circuit . 

In the data interpolating device, a plurality of 
over- sampling circuits may be cascade connected, when one set 
of over-sampling circuit consists of the plural stages of delay 
circuits, the multiplication/addition circuit and the 
averaging circuit . 

Also, the invention provides a data interpolating device 
comprising four stages of delay circuits for successively 
delaying discrete data sequentially inputted, and a 
multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the four 
stages of delay circuits according to the value of a digital 
basic function and thereby determines interpolation data for 
the discrete data sequentially inputted. 

In the data interpolating device, a plurality of 
over- sampling circuits may be cascade connected, when one set 
of over- sampling circuit consists of the four stages of delay 
circuits and the multiplication/addition circuit. 

Also, the invention provides a data interpolating device 
comprising five stages of delay circuits for successively 
delaying discrete data sequentially inputted, and a 
multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the five 
stages of delay circuits according to the value of a digital 
basic function and thereby determines interpolation data for 
the discrete data sequentially inputted. 
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In the data interpolating device, a plurality of 
over- sampling circuits may be cascade connected, when one set 
of over- sampling circuit consists of the five stages of delay 
circuits and the multiplication/addition circuit. 

Also, the invention provides a data interpolating device 
comprising data acquisition means for acquiring discrete data 
at a sampling point of notice and discrete data at neighboring 
sampling points around the sampling point of notice when the 
discrete data are sequentially inputted, and interpolation 
means that performs weighted addition of discrete data acquired 
by the data acquisition means according to the value of a digital 
basic function and thereby determines interpolation data for 
the discrete data at the sampling point of notice successively . 

In the data interpolating device, the operation of the 
data acquisition means and the interpolation means may be 
repeated multiple times by setting the output data of the 
interpolation means as the input of the data acquisition means . 

The data interpolating device may further comprise 
averaging means for performing an averaging operation of 
adjacent interpolation data for the interpolation data 
determined by the interpolation means. 

The invention provides a data interpolating method 
comprising performing weighted addition of discrete data at 
a sampling point of notice and discrete data at surrounding 
sampling points around the sampling point of notice according 
to the value of a digital basic function and thereby determining 
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interpolation data for the discrete data at the sampling point 
of notice, when the discrete data are sequentially inputted. 

The data interpolating method may further comprise 
performing an averaging operation of adjacent interpolation 
data for the interpolation data determined by performing 
weighted addition of the discrete data according to the value 
of the digital basic function. 

The data interpolating method may further comprise 
performing weighted addition of interpolation data at a 
sampling point of notice and interpolation data at surrounding 
sampling points around the sampling point of notice according 
to the value of the digital basic function, when the 
interpolation data are determined by performing weighted 
addition of discrete data according to the value of the digital 
basic function , and thereby further determining interpolation 
data for the interpolation data at the sampling point of notice . 

The invention provides a data interpolating method, 
comprising performing weighted addition of discrete data at 
a sampling point of notice and discrete data at its adjacent 
sampling points according to the value of a digital basic 
function when the discrete data are sequentially inputted and 
thereby successively determining interpolation data at two 
sampling points from the discrete data at the one sampling 
point . 

Also, the invention provides a data interpolating method 
comprising performing weighted addition of discrete data at 
two adjacent sampling points of notice and discrete data at 
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its adjacent sampling points of the two sampling points 
according to the value of a digital basic function when the 
discrete data are sequentially inputted and thereby 
successively determining interpolation data interpolating 
between the two sampling points. 

This invention provides a sampling function generating 
device comprising plural stages of delay circuits for 
successively delaying singly inputted discrete data, a 
multiplication/addition circuit that performs weighted 
addition of data outputted from the output stages of the plural 
stages of delay circuits according to the value of a digital 
basic function and thereby determines interpolation data for 
the singly inputted discrete data, and a plurality of 
over- sampling circuits that are cascade connected when one 
set of over-sampling circuit consists of the plural stages 
of delay circuits and the multiplication/addition circuit. 

The invention provides a data interpolation program for 
enabling a computer to operate as each means according to claim 
17, or to perform a processing procedure of the data 
interpolation method according to any one of claims 20, 24 
and 26. 

Also , the invention provides a computer readable recorded 
medium storing a program for enabling a computer to operate 
as each means according to claim 17 , or to perform a processing 
procedure of the data interpolation method according to any 
one of claims 20, 24 and 26. 



[Brief Description of the Drawings] 

FIG. 1 is an explanatory view of a sine function; 

FIG. 2 is an explanatory view of data interpolation using 
a sampling function as shown in FIG. 1; 

FIG. 3 is a characteristic graph showing a conventional 
sampling function generated from a digital basic function as 
shown in FIG. 15; 

FIG. 4 is a diagram showing a configuration example of 
a data interpolating device according to a first embodiment; 

FIG. 5 is a diagram for explaining a principle of data 
interpolation according to the first embodiment; 

FIG. 6 is a characteristic graph showing the over- sampled 
result obtained when unit pulse data is input into the data 
interpolating device as shown in FIG. 4; 

FIG. 7 is a characteristic graph showing the comparison 
of a waveform function as indicated by #4 in FIG. 6 and a waveform 
function as shown in FIG. 3; 

FIG. 8 is a graph, partially enlarged, of a sampling 
function as shown in FIG. 7; 

FIG. 9 is a graph representing the frequency 
characteristic of the sampling function as shown in FIG. 7; 

FIG . 10 is a diagram showing another configuration example 
of the data interpolating device according to the first 
embodiment ; 

FIG. 11 is a diagram showing another configuration example 
of the data interpolating device according to the first 
embodiment; 
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FIG . 12 is a diagram showing another configuration example 
of the data interpolating device according to the first 
embodiment; 

FIG. 13 is a diagram for explaining a principle of data 
interpolation according to a second embodiment; 

FIG. 14 is a diagram showing a configuration example of 
the data interpolating device according to the second 
embodiment ; and 

FIG. 15 is a graph showing a digital basic function for 
use with this embodiment. 

[Detailed Description of the Preferred Embodiments] 
(First embodiment) 

A first embodiment of the present invention will be 
described below with reference to the accompanying drawings . 

FIG. 4 is a diagram showing a configuration example of 
a data interpolating device according to the first embodiment . 
FIG. 5 is a diagram for explaining a principle of data 
interpolation according to the first embodiment. Referring 
firstly to FIG. 5, the principle of data interpolation 
according to this embodiment will be described below. 

In a data structure as shown in FIG. 5, A, B, C, ... along 
the transverse axis indicate discrete data at each sampling 
point sequentially inputted at every clock, in which its data 
value becomes a coefficient of a digital basic function . Also , 
a, b, c, ... along the longitudinal axis indicate a central 



position of discrete data processed by the digital basic 
function. 

The digital basic function as used herein is a base for 
a sampling function useful in making data interpolation by 
over- sampling, as shown in FIG. 13. This digital basic 
function is produced by changing data value such as -1, 1, 
8 # 8, 1 and -1 at every clock. 

As shown in FIG. 5, the digital basic function for discrete 
data A is given by shifting each function value (-1, 1, 8, 
8,1, -1) by every one clock from the start point that is the 
first clock position. Also, the digital basic function for 
the next discrete data B is given by shifting each function 
value by every one clock from the start point that is a clock 
position two clocks after the first clock position . Similarly, 
the digital basic function for discrete data C, D, E, F, ... 
sequentially inputted is given by shifting each function value 
by every one clock from the start point that is a clock position 
more two clocks delayed . 

An instance where two interpolation data Bl, B2 are 
produced from one discrete data B on the basis of this data 
structure. Herein, assuming that two lines of data across 
the central position b of the discrete data B processed by 
the digital basic function are bl and b2 , the convolution 
operation of the data bl and b2 is performed. In this case, 
two adjacent data values bl and b2 are represented such as 
bl=A+8B-C ... (l) 

b2=-A+8B+C ... (2) 
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From the above expressions (1) and (2), 
bl+b2=16B ... (3) 

results, and is further transformed such as 
B=(bl/8+b2/8)/2 ... (4) 

From the above expression (4) , it will be found that the 
discrete data B corresponds to an intermediate value between 
two data bl/8 and b2/8. Conversely, one discrete data B is 
replaced with two interpolation data Bl and B2, which are 
represented such as 

Bl=(A+8B-C)/8 ... (5) 

B2=(-A+8B+C)/8 ... (6) 

Since the relation of expression (4) for the discrete 
data B similarly holds for other discrete data C, D, E, F, 
C=(cl/8+c2/8)/2 
D=(dl/8+d2/8)/2 
E=(el/8+e2/8)/2 
F=(f l/8+f2/8)/2 

are obtained. 

Accordingly, one discrete data C is replaced with two 
interpolation data CI and C2, one discrete data D is replaced 
with two interpolation data Dl and D2, one discrete data E 
is replaced with two interpolation data El and E2, and one 
discrete data F is replaced with two interpolation data Fl 
and F2, as shown in the following expressions. 

C -> Cl=(B+8C-D)/8, C2=(-B+8C+D)/8 

D -* Dl=(C+8D-E)/8, D2=(-C+8D+E)/8 
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E -» El=(D+8E-F)/8 f E2= ( -D+8E+F ) /8 
F -» Fl=(E+8F-G)/8, F2= ( -E+8F+G ) /8 

• ••»». 

As above described, when the interpolation value is 
obtained for the discrete data (e.g. , B) at a certain sampling 
point, the value of the discrete data at the sampling point 
multiplied by 8 and the values of the discrete data at previous 
and next sampling points multiplied by +1 and -1 are added, 
and the added value is divided by 8 to obtain a first interpolated 
value (e.g., Bl). Also, the value of the discrete data of 
object multiplied by 8 and the values of the discrete data 
at the previous and next sampling points multiplied by -1 and 
+ 1 are added, and the added value is divided by 8 to obtain 
a second interpolation value (e.g. , B2 ) . By performing this 
operation for each sampling point, the original data is 
over- sampled twice. 

Also , if the same interpolating process is performed based 
on the interpolation data over- sampled twice, the original 
data is over- sampled four times. Moreover, if the same 
interpolating process is performed based on the interpolation 
data over- sampled four times , the original data is over- sampled 
eight times. Theoretically, if the same interpolating 
process is repeated n times , the original data is over- sampled 
2 n times. 

Referring to FIG. 4, the configuration of the data 
interpolating device for implementing the interpolating 
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operation according to the first embodiment will be described 
below. 

The data interpolating device as shown in FIG. 4 is 
configured such that three stages of two-fold over-sampling 
circuits are cascade connected, in which the circuit 
configurations at the first to third stages are distinguished 
by subscript (-1 to -3). Also, the circuit configuration 
having the same main symbol with different subscript has the 
same function. 

The first stage over - sampling circuit operates in 
accordance with the clock CK of a reference frequency (e.g. , 
44.1 KHz). Three stages of D-type flip-flops l_i, 2_ x and 3-! 
delay discrete data inputted in sequence ( for example , 16 bits ) 
every one clock CK of a reference frequency sequencially . 
These D-type flip-flops l_i, 2. 1# and 3_! correspond to three 
stages of delay circuits according to the invention. 

Data taken out from an output tap of the first stage D-type 
flip-flop l_ x forms the first term as represented by the 
expressions (5) and (6), data taken out from an output tap 
of the second stage D-type flip-flop 2_ x forms the second term, 
and data taken out from an output tap of the third stage D-type 
flip-flop 3>i forms the third term, for example. 

Data taken out from the output tap of the first stage 
D-type flip-flop l-i is entered via a -1 time multiplication 
unit 4_! (corresponding to the first multiplication unit as 
defined in claim 3) into one input terminal of an AND gate 
5_i, and entered into one input terminal of an AND gate 6-i 



- 17 - 



without interposition of the multiplication unit 
(corresponding to the +1 time). A reverse clock CK passing 
through an inverter 7-! is input at the other input end of 
the AND gate 5_i. Also, clock CK is input at the other input 
end of an AND gate 6_i. 

The output data of the two AND gates 5_i and 6_ x are output 
via an OR gate 8_i. Thereby, the discrete data multiplied by 
+ 1 is output from the OR gate 8_i during the period for which 
the clock CK is at M H" . Also, the discrete data multiplied 
by -1 is output from the OR gate 8-i during the period for 
which the clock CK is at "L" . That is, the first term of the 
expression (5) is obtained during the period for which the 
clock CK is at "H" , and the first term of the expression (6) 
is obtained during the period for which the clock CK is at 
"L n . A first switching circuit as defined in claim 3 is 
composed of the two AND gates 5-! and 6. 1# the inverter 7_i and 
the OR gate 8.!. 

Data taken out from the output tap of the second stage 
D-type flip-flop 2- X is entered via a +8 times multiplication 
unit 9_i (corresponding to the second multiplication unit as 
defined in claim 3). Since the second term is not changed 
in ± sign in both the expressions (5) and- (6), there is no 
need for a sign switching circuit with the clock CK as in the 
first term. 

Also, data taken out from the output tap of the third 
stage D-type flip-flop 3_i is entered via a -1 time 
multiplication unit 10.! (corresponding to the third 
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multiplication unit as defined in claim 3) into one input 
terminal of an AND gate 11. i. and entered into one input terminal 
of an AND gate 12. x without interposition of the multiplication 
unit (corresponding to the +1 time) . A clock CK is input into 
the other input end of the AND gate 11.x . Also . a reverse clock 
CK passing through an inverter 13. x is input at the other input 
end of the AND gate 12.! . 

The output data of the two AND gates ll-i and 12_i are 
output via an OR gate 14_i. Thereby, the discrete data 
multiplied by -1 is output from the OR gate 14. x during the 
period for which the clock CK is at "H". Also, the discrete 
data multiplied by +1 is output from the OR gate 14. x during 
the period for which the clock CK is at "L". That is, the 
third term of the expression ( 5 ) is obtained during the period 
for which the clock CK is at "H", and the third term of the 
expression (6) is obtained during the period for which the 
clock CK is at "L" . A second switching circuit as defined 
in claim 3 is composed of the two AND gates 11.! and 12. i. the 
inverter 13. x and the OR gate 14- x . 

The output data of the OR gate 8_i, the output data of 
the 8 times multiplication unit 9_i and the output data of 
the OR gate 14.! are added by two adders 15. 1 and 16. i 
(corresponding to the adder as defined in claim 3) . Thereby, 
the operation result of the expression ( 5 ) is output from adder 
16.1 during the period for which the clock CK is at "H" , and 
the operation result of the expression ( 6 ) is output from adder 
16-i during the period for which the clock CK is at "L". 
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If the discrete data A, B, C, ... as shown in FIG. 5 is 
sequentially input into the above circuit, the interpolation 
data over-sampled twice such as A, Bl, B2, CI, C2, ... are 
output from the adder 16. i. Though the data value is wholly 
multiplied by 1/8 in the expressions (5) and (6), no circuit 
therefor is provided in the first stage over- sampling circuit . 
In the circuit of FIG. 4, data at three stages are multiplied 
by 1/512 ( 1/8 3 ) in a multiplication unit 20 at the final stage. 

In the data interpolating device as shown in FIG. 4, when 
the interpolation data Bl and B2 are produced from one discrete 
data B, for example, the second terms of the expressions (5) 
and (6) common in both the interpolation data share one 8 times 
multiplication unit 9-i to form the data. In this manner, the 
circuit to be shared in producing two interpolation data is 
shared if possible, simplifying the entire circuit 
configuration . 

An averaging circuit composed of the D-type flip-flop 
17_i, an adder 18-! , and a 1/2 times multiplication unit 19-! 
is provided at the output stage of the adder 16-! . The D-type 
flip-flop 17-i delays the output data of the adder 16_i by one 
clock in accordance with the clock 2CK that is twice the 
reference frequency. The adder 18-i adds the output data of 
the adder 16_i and the output data of the D-type flip-flop 17-!. 
The 1/2 times multiplication unit 19_i multiplies the output 
data of adder 18-i by 1/2. 

The reason of providing the averaging circuit is to restore 
a deviation in the relative position of clock caused by the 
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interpolation process. For example, when the interpolation 
data Bl, B2, CI, C2 , . .. are output from the adder 16. i# the 
averaging circuit performs the operations (Bl+B2)/2, 
(B2+Cl)/2, (Cl+C2)/2, ... successively. 

The execution result of this averaging operation is as 
follows . 

(Bl+B2)/2={{(A+8B-C)/8}+{(-A+8B+C)/8}}/2=B 
(B2+Cl)/2={{ (-A+8B+C)/8}+{ ( B+8C-D ) /8 } >/2 

=(-A+9B+9C-D)/16 ... (7) 

(Cl+C2)/2={{(B+8C-D)/8}+{(-B+8C+D)/8}}/2=C 

As a result, the deviation in the relative position of clock 
is restored to the original correct position. Thereby, it 
is possible to make better the frequency characteristics of 
the obtained sampling function as will be described later with 
reference to FIG. 9. 

The output data of 1/2 times multiplication unit 19-! 
composing the averaging circuit is input into the D-type 
flip-flop 1. 2 of the second stage over-sampling circuit, so 
that the two-fold over-sampling at the second stage is 
performed. The over- sampling circuit at the second stage is 
the same as that at the first stage, except that it operates 
in accordance with the clock 2CK that is twice the reference 
frequency. 

The output data of 1/2 times multiplication unit 19. i 
provided at the final stage of the second stage over-sampling 
circuit is input into the D-type flip-flop 1_ 3 of the third 
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stage over-sampling circuit, so that the two-fold 
over- sampling at the third stage is performed. The 
over- sampling circuit at the third stage is the same as that 
at the first stage, except that it operates in accordance with 
the clock 4CK that is four times the reference frequency. 

FIG. 6 is a characteristic graph showing the over- sampled 
result obtained when unit pulse data is input into the data 
interpolating device as shown in FIG. 4. In FIG. 6, #1 input 
data of unit pulse, #2 two-f old over- sampled data, #3 four-fold 
over-sampled data, and #4 eight-fold over-sampled data are 
shown . A waveform function of over- sampled data as indicated 
by #4 has the almost same waveform as the sampling function 
obtained by the convolution operation as described in Japanese 
Patent Application No. 11-173245, as shown in FIG. 3. 

For reference, the waveform function as indicated by #4 
in FIG. 6 and the waveform function as shown in FIG. 3 are 
shown on the same plane for comparison in FIG. 7. In FIG. 
7, graph I of solid line is the waveform function as indicated 
by #4 in FIG. 6 and graph II of dotted line is the waveform 
function as shown in FIG. 3. As will be apparent from these 
graphs, there is almost no difference in both the waveforms, 
and these graphs represent the sampling function of finite 
supports converging to zero at finite sampling positions . And 
the function as indicated by #4 in FIG. 6 is dif f erentiable 
once or more times over the whole range. 

Accordingly, in the case where the data interpolating 
process of this embodiment is applied, superposition is made 
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based on each discrete data, employing the sampling function 
as indicated by #4 in FIG. 6 as a result, whereby it is possible 
to interpolate more smoothly the value between each discrete 
data, employing the first or higher order different iable 
function. Thereby, the low pass filter is dispensed with, 
and the degradation in the phase characteristic is prevented. 

Also, in the case where the interpolation value for one 
discrete data is obtained, it is only necessary to take into 
consideration the values for a limited number of discrete data 
( three as represented in the expressions ( 5 ) and ( 6 ) ) , whereby 
it is possible to reduce the processing amount significantly. 
Though each discrete data outside the range of finite supports 
should be naturally considered, but are not ignored in 
consideration of the processing amount or precision. 
Theoretically, there is no need for considering them, whereby 
it is possible to prevent a truncation error from occurring. 

FIG. 8 is a graph, partially enlarged, of the sampling 
function as shown in FIG. 7. As shown in FIG. 8, the sampling 
function 11 acquired from Japanese Patent Application No. 
11-173245 is compulsorily converged to zero at sampling point 
t=7. However, the sampling function I in this embodiment is 
not compulsorily truncated. Accordingly, according to this 
embodiment , it is possible to interpolate between each discrete 
data more smoothly. 

FIG. 9 is a graph representing the frequency 
characteristic of the sampling function as shown in FIG. 7. 
In FIG. 9, graph I of solid line represents the frequency 
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characteristics of the sampling function acquired in this 
embodiment, and graph II of dotted line represents the 
frequency characteristics of the sampling function acquired 
from Japanese Patent Application No. 11-173245. As will be 
seen from FIG. 9, the extremely excellent frequency 
characteristics are obtained in both graphs . Thereby, a notch 
filter is placed where a so-called reflected noise appears, 
whereby the reflected noise is effectively suppressed. 

As described above in detail, according to the first 
embodiment, the sampling function of finite supports 
dif f erentiable once or more times over the whole range is 
obtained only by very simple arithmetical operations as 
represented in the expressions (5) and (6), and data 
interpolation based on the sampling function is made by simple 
arithmetical operations . In the expressions ( 5 ) and ( 6 ) , when 
interpolation data is obtained for discrete data at a sampling 
point of notice, it is necessary to take into consideration 
only the discrete data at its adjacent sampling points. 
Accordingly, according to this embodiment, it is possible to 
not only obtain excellent interpolation data without phase 
distortion or truncation error of the low pass filter, but 
also shorten the operation time and simplify the operation 
circuit . 

The data interpolating circuit as shown in FIG. 4 is one 
example of the circuit configuration for implementing the data 
interpolating method of the invention, but is not limited to 
this example. 
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For example, a first multiplication/addition circuit for 
adding the data output from the output taps of three D-type 
flip-flops and multiplied by -1, 8 and +1, and a second 
multiplication/addition circuit for adding the data output 
from the output taps of three D-type flip-flops and multiplied 
by +1, 8 and -1 may be provided apartly. 

FIG. 10 is a diagram showing a circuit configuration 
example in this case. Though a two-fold over- sampling circuit 
is only shown, the n same circuits may be cascade connected 
to make 2 n -fold over- sampling in the same manner as in FIG. 
4. The averaging circuit is not shown herein. 

In FIG. 10, three stages of D-type flip-flops 21, 22 and 
23 delay the sequentially inputted discrete data by every one 
clock CK successively. A first multiplication/addition 
circuit comprises a first multiplication unit 24 for 
multiplying the output data of a first stage D-type flip-flop 
21 by -1, a second multiplication unit 25 for multiplying the 
output data of a second stage D-type flip-flop 22 by 8, and 
the adders 26 and 27 for adding the output data of the first 
multiplication unit 24, the output data of the second 
multiplication unit 25, and the output data of a third stage 
D-type flip-flop 23. 

Also , a secondmultiplication/addition circuit comprises 
a third multiplication unit 28 for multiplying the output data 
of the second stage D-type flip-flop 22 by 8, a fourth 
multiplication unit 29 for multiplying the output data of the 
third stage D-type flip-flop 23 by -l f and the adders 30 and 
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31 for adding the output data of the third multiplication unit 
28, the output data of the fourth multiplication unit 29, and 
the output data of the first stage D-type flip-flop 21. 

Data outputted from the first multiplication/addition 
circuit (adder 27) is entered into one input terminal of an 
AND gate 32. A reverse clock CK passing through an inverter 
33 is input at the other input end of the AND gate 32. Also, 
data outputted from the second multiplication/addition 
circuit (adder 31) is entered into one input terminal of an 
AND gate 34. A clock CK is input at the other input end of 
the AND gate 34. 

The output data of the two AND gates 32 and 34 are output 
via an OR gate 35 . Thereby, the operation result as indicated 
in the expression (5) is output from the OR gate 35 during 
the period for which the clock CK is at "H" . Also , the operation 
result as indicated in the expression (6) is output from the 
OR gate 35 during the period for which the clock CK is at "L" . 

Also, though two eight-fold multiplication units are 
employed in FIG. 10, a single multiplication unit may be 
employed as shown in FIG. 11 . Though a two -fold over- sampling 
circuit is only shown in FIG. 11, the n same circuits may be 
cascade connected to make 2 n -fold over- sampling . The 
averaging circuit is not shown herein. 

The data interpolating device as shown in FIG . 11 comprises 
three stages of D-type flip-flops 41, 42 and 43 for delaying 
the sequentially inputted discrete data by every one clock 
CK successively. Also, it comprises a first multiplication 
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unit 44 for multiplying the output data of a first stage D-type 
flip-flop 41 by -1, a second multiplication unit 45 for 
multiplying the output data of a second stage D-type flip-flop 
42 by 8, and a third multiplication unit 46 for multiplying 
the output data of a third stage D-type flip-flop 43 by -1. 

Further, it comprises the first adders 47 and 48 for adding 
the output data of the first multiplication unit 44 , the output 
data of the second multiplication unit 45, and the output data 
of the third stage D-type flip-flop 43 and the second adders 
49 and 50 for adding the output data of the secondmultiplication 
unit 45, the output data of the third multiplication unit 46, 
and the output data of the first stage D-type flip-flop 41. 

Data outputted from an adder 48 is entered into one input 
terminal of an AND gate 51 . A reverse clock CK passing through 
an inverter 52 is input at the other input end of the AND gate 
51, Also, data outputted from an adder 50 is entered into 
one input terminal of an AND gate 53. A clock CK is input 
at the other input end of the AND gate 53. 

The output data of the two AND gates 51 and 53 are output 
via an OR gate 54 . Thereby, the operation result as indicated 
in the expression (5) is output from the OR gate 54 during 
the period for which the clock CK is at "H" . Also, theoperation 
result as indicated in the expression (6) is output from the 
OR gate 54 during the period for which the clock CK is at n L n . 

In the above embodiment, data interpolation is performed 
by making the arithmetical operation as shown in the 
expressions (5) and (6), and a deviation in the relative 
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position of clock is compensated by making the averaging 
operation of obtained interpolation data. On the contrary, 
the averaging operation may be omitted by directly making the 
arithmetical operation as shown in the expression (7). 

When the arithmetical operation as shown in the expression 
(7) is directly made, the weighted addition of discrete data 
at two adjacent sampling points of notice and discrete data 
at directly adjacent sampling points is performed according 
to the value of a digital basic function when the discrete 
data are sequentially inputted and thereby interpolation data 
interpolating the discrete data between two sampling points 
is successively determined. For example, when two discrete 
data of notice are B and C, interpolation data (B+C)/2 between 
them is obtained as in the following expression. 
(B+C)/2=(-A+9B+9C-D)/16 . (8) 

FIG. 12 is a diagram showing a circuit configuration 
example in this case . Though a two -fold over- sampling circuit 
is only shown here , the n same circuits may be cascade connected 
to make 2 n -fold over- sampling. 

The data interpolating device as shown in FIG. 12 comprises 
four stages of D-type flip-flops 61, 62 63 and 64 for delaying 
the sequentially inputted discrete data by every one clock 
CK successively. The four stages of D-type flip-flops 61 to 
64 correspond to the four stages of delay circuits according 
to the invention. 

Also, it comprises a first multiplication unit 65 for 
multiplying the output data of a first stage D-type flip-flop 
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61 by -1, a second multiplication unit 6 6 for multiplying the 
output data of a second stage D-type flip-flop 62 by 9, a third 
multiplication unit 67 for multiplying the output data of a 
third stage D-type flip-flop 63 by 9, and a fourth 
multiplication unit 68 for multiplying the output data of a 
fourth stage D-type flip-flop 64 by -1 . Moreover , it comprises 
the adders 69, 70 and 71 for adding all the output data of 
the first to fourth multiplication units 65 to 68. 

Data output ted from an adder 71 is entered into one input 
terminal of an AND gate 72 . A reverse clock CK passing through 
an inverter 73 is input at the other input end of the AND gate 
72 . Also , the discrete data entered into the first stage D-type 
flip-flop 61 is input via a delay circuit 74 for delaying data 
by an amount corresponding to a delay of each circuit block 
61 to 71 at one input terminal of the AND gate 75. A clock 
CK is input at the other input end of the AND gate 75. 

The output data of the two AND gates 72 and 75 are output 
via an OR gate 76 . Thereby, the input discrete data is directly 
output from the OR gate 76 during the period for which the 
clock CK is at "H" . Also, the interpolation data as indicated 
in the expression (8) is output from the OR gate 76 during 
the period for which the clock CK is at "L". 

FIG . 12 is one circuit example for performing the operation 
of the expression (8), but the circuit for performing the 
operation of the expression (8) is not limited to this example. 
For example, the following configuration may be employed 
instead of the multiplication/addition circuits 65 to 71 as 
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shown in FIG. 12. That is, the multiplication/addition 
circuit comprises a first adder for adding the output data 
of the first stage D-type flip-flop 61 and the output data 
of the fourth stage D-type flip-flop 64, a second adder for 
adding the output data of the second D-type flip-flop 62 and 
the output data of the third stage D-type flip-flop 63 , a first 
multiplication unit for multiplying the output data of the 
first adder by -1, a second multiplication unit for multiplying 
the output data of the second adder by 9, and a third adder 
for adding the output data of the first multiplication unit 
and the output data of the second multiplication unit. In 
this manner, the number of multiplication units is reduced. 
(Second embodiment) 

A second embodiment of the invention will be described 
below. In the second embodiment, a function in which the 
digital basic function of FIG. 15 shifted by every one clock 
is added (averaged) is employed. 

FIG. 13 is a diagram for explaining a principle of data 
interpolation according to the second embodiment . In a data 
structure as shown in FIG. 13, A, B, C, ... along the transverse 
axis indicate discrete data at each sampling point sequentially 
inputted at every clock, in which its data value becomes a 
coefficient of a digital basic function in this embodiment. 
Also, a, b, c, ... along the longitudinal axis indicate a central 
position of discrete data processed by the digital basic 
function. 
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As shown in FIG, 13, the digital basic function for 
discrete data A is given by shifting each function value (-1, 
0, 9, 16, 9, 0, -1) by every one clock from the start point 
that is two clocks ahead of the first clock position. Also, 
the digital basic function for the next discrete data B is 
given by shifting each function value by every one clock from 
the start point that is the first clock position. Similarly, 
the digital basic function for discrete data C, D, E f F, ... 
sequentially inputted is given by shifting each function value 
by every one clock from the start point that is a clock position 
more two clocks delayed. 

With this data structure, the central position c of the 
digital basic function processed according to the value of 
discrete data C lies on one data position, unlike the first 
embodiment. Thus, in the second embodiment, the convolution 
operation is performed employing data at the clock position 
c and data at the clock positions before and after the clock 
position c to determine two interpolation data CI and C2 from 
one thinned data C. 

Supposing three lines of data located at the clock position 
c and before and after it are c- , c and c+ , they are represented 
such as 

c-=-A+9B+9C-D ... (9) 
c=16C ... (10) 

C+=-B+9C+9D-E ... (11) 

From the above expressions (9) and (10), 
(c-+c)=(-A+9B+25C-D) ... (12) 
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is obtained. And from the above expressions (10) and (11), 

(c+c+)=(-B+25C+9D-E) . . . (13) 
is obtained. 

Taking into consideration an amount of adjusting the 
amplitude based on the expressions ( 12 ) and ( 13 ) , one discrete 
data is replaced with two interpolation data CI and C2 , which 
are represented such as 

Cl=(-A+9B+25C-D)/32 . . . (14) 

C2=(-B+25C+9D-E)/32 . . . (15) 

Also, the relation as represented by the expressions (14) 
and (15) for the discrete data C holds for other discrete data 
D, E # F, ... Accordingly, one discrete data D is replaced with 
two interpolation data Dl and D2 , one discrete data E is replaced 
with two interpolation data El and E2, and one discrete data 
F is replaced with two interpolation data Fl and F2, as shown 
in the following expression. 

D -> Dl=(-B+9C+25D-E)/32, D2= ( -C+25D+9E-F ) /32 
E El=(-C+9D+25E-F)/32, E2= ( -D+25E+9F-G) /32 
F Fl=(-D+9E+25F-G)/32, F2= ( -E+25F+9G-H) /32 
... ... 

As above described, when the interpolation value is 
obtained for the discrete data (e.g. , C) at a certain sampling 
point, the value of its discrete data multiplied by 25, the 
value of the discrete data at the previous sampling point 
multiplied by 9, the value of the discrete data ahead of two 
sampling points multiplied by - 1 , and the value of the discrete 
data at next sampling point multiplied by -1 are added, and 
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the added value is divided by 8 to obtain a first interpolated 
value (e.g., CI). 

Also, the value of the discrete data of object multiplied 
by 25, the value of the discrete data at the previous sampling 
point multiplied by -1, the value of the discrete data at the 
next sampling point multiplied by 9, and the value of the 
discrete data after two sampling points multiplied by -1 are 
added, and the added value is divided by 8 to obtain a second 
interpolation value (e.g., C2). 

By performing this operation for each sampling point, 
the original data is over- sampled twice. 

Also , if the same interpolating process is performed based 
on the interpolation data over- sampled twice, the original 
data is over-sampled four times. Moreover, if the same 
interpolating process is performed based on the interpolation 
data over- sampled four times , the original data is over- sampled 
eight times. Theoretically, if the same interpolating 
process is repeated n times , the original data is over- sampled 
2 n times. 

FIG. 14 is a diagram showing a configuration example of 
the data interpolating device according to the second 
embodiment. The data interpolating device as shown in FIG. 
14 is configured such that three stages of two-fold 
over- sampling circuits are cascade connected, in which the 
circuit configurations at the first to third stages are 
distinguished by subscript (-1 to -3). Also, the circuit 
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configuration having the same main symbol with different 
subscript has the same function. 

The first stage over- sampling circuit operates in 
accordance with the clock CK of a reference frequency (e.g. , 
44.1 KHz). Five stages of D-type flip-flops 81. 1# 82. 1# 83. lt 
84.! and 85. i delay the sequentially inputted discrete data 
(e.g., 16 bits ) by every one clock CK of the reference frequency 
successively. These D-type flip-flops 81-i to 85_i correspond 
to five stages of delay circuits as defined in claim 13 of 
the invention. 

Data taken out from the output tap of the first stage 
D-type flip-flop 81_! is entered into a -1 time multiplication 
unit 86. x (corresponding to the first multiplication unit as 
defined in claim 14). Data taken out from the output tap of 
the second stage D-type flip-flop 82-! is entered into a +9 
times multiplication unit 87. i (corresponding to the second 
multiplication unit as defined in claim 14) and a -1 time 
multiplication unit 93_i (corresponding to the fifth 
multiplication unit as defined in claim 14). 

Data taken out from the output tap of the third stage 
D-type flip-flop 83. x is entered into a 25 times multiplication 
unit 88-i (corresponding to the third multiplication unit as 
defined in claim 14). Data taken out from the output tap of 
the fourth stage D-type flip-flop 84.! is entered into a -1 
time multiplication unit 89. i (corresponding to the fourth 
multiplication unit as defined in claim 14) and a +9 times 
multiplication unit 94_i (corresponding to the sixth 
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multiplication unit as defined in claim 14) . Data taken out 
from the output tap of the fifth stage D-type flip-flop 85_i 
is entered into a -1 time multiplication unit 95. i 
(corresponding to the seventh multiplication unit as defined 
in claim 14 ) . 

The output data of the -1 time multiplication unit 86. i, 
the 9 times multiplication unit 87. lt the 25 times 
multiplication unit 88. x and the -1 time multiplication unit 
89_i are all added by three adders 90-x , 91_i and 92_i . The output 
data of the -1 time multiplication unit 93. i# the 25 times 
multiplication unit 88. lt the 9 times multiplication unit 94_ x 
and the -1 time multiplication unit 95_i are all added by three 
adders 96-i, 97. i and 98_i. 

Data output from the adder 92_i is input into one input 
terminal of an AND gate 99. x . A clock CK is input at the other 
input end of an AND gate 9 9-i. Also, data output from the adder 
98.i is input into one input terminal of an AND gate 100. i. 
A reverse clock CK passing through an inverter 101. i is input 
at the other input end of the AND gate 100. i. 

The output data of the two AND gates 99-i and lOO-i are 
output via an OR gate 102_i. Thereby, the operation result 
of the expression (14) is output from the OR gate 102. i during 
the period for which the clock CK is at "H" , and the operation 
result of the expression (15) is output from the OR gate 102. x 
during the period for which the clock CK is at n L" . 

If the discrete data A, B, C, ... as shown in FIG. 13 
is sequentially input into the above circuit , the interpolation 
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data over- sampled twice such as A, B , CI , C2 , Dl , D2 , ... are 
output from the OR gate 102^. Though the data value is wholly 
multiplied by 1/32 in the expressions ( 14 ) and ( 15 ) , no circuit 
therefor is provided in the first stage over-sampling circuit . 
In the circuit of FIG. 14, data at three stages are multiplied 
by 1/32768 (1/32 3 ) in a multiplication unit 104 at the final 
stage . 

The output data of the OR gate 102-i is input into the 
D-type flip-flop 81. 2 of the second stage over-sampling circuit , 
so that the two-fold over-sampling at the second stage is 
performed. The over-sampling circuit at the second stage is 
the same as that at the first stage, except that it operates 
in accordance with the clock 2CK that is twice the reference 
frequency. 

The output data of the OR gate 102. 2 provided at the final 
stage of the second stage over- sampling circuit is input into 
the D-type flip-flop 81. 3 of the third stage over-sampling 
circuit, so that the two-fold over-sampling at the third stage 
is performed. The over- sampling circuit at the third stage 
is the same as that at the first stage, except that it operates 
in accordance with the clock 4CK that is four times the reference 
frequency. 

The output data of the OR gate 102. 3 provided at the final 
stage of the three stage over- sampling circuit is held in the 
D-type flip-flop 103 in accordance with the clock 8CK of 
eight-fold frequency, and then output via the 1/512 times 
multiplication unit 104. 
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When the data interpolating device is configured as shown 
in FIG. 14, the sampling function obtained in inputting the 
unit pulse data is the almost same waveform as indicated by 
#4 in FIG. 6. Accordingly, in the case where the data 
interpolating process of the second embodiment is applied,, 
it is possible to interpolate more smoothly the value between 
each discrete data, employing the first or higher order 
different iable function. Thereby, the low pass filter is 
dispensed with , and the degradation in the phase characteristic 
is prevented. 

Also, in the case where the interpolation value for one 
discrete data is obtained, it is only necessary to take into 
consideration the values for a limited number of discrete data 
( four as represented in the expressions ( 14 ) and ( 15 ) ) , whereby 
it is possible to reduce the processing amount significantly. 
Though each discrete data outside the range of finite supports 
should be naturally considered, but are not ignored in 
consideration of the processing amount or precision. 
Theoretically, there is no need for considering them, whereby 
it is possible to prevent a truncation error from occurring. 

Also, in the second embodiment, interpolation data is 
obtained employing a function in which the digital basic 
function of FIG. 15 shifted by every one clock is added 
(averaged) as shown in FIG. 13 , whereby the averaging operation 
is already performed in this function. Accordingly, there 
is no need for performing the averaging operation after the 
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interpolating operation in accordance with the expressions 
(14) and (15), whereby the averaging circuit is omitted. 

In the first and second embodiments, the digital basic 
function has the data value changed at every clock such as 
■l.i 8 * 8 * land-1. Though these numerical values are optimal, 
the digital basic function may have other numerical values 
(e.g., weight at both ends is not -1 but 1 or 0, weight at 
the center is other than 8) . Even in this case, the excellent 
interpolation data interpolating smoothly between each 
discrete data is obtained through the simple operation as 
compared with when the sine function is conventionally 
employed . 

The method of data interpolation according to the first 
and second embodiments as described above may be implemented 
by any of the hardware, DSP and software configuration. For 
example, when this method is implemented by software, the data 
interpolating device of this embodiment is actually composed 
of CPU, MPU, RAM and ROM in the computer to execute a program 
stored in the RAM or ROM. 

Accordingly, the method is implemented by storing the 
program for enabling the computer to perform the functions 
of the above embodiments in the recording medium such as CD-ROM, 
and read the program when executed. The recording medium for 
recording the program may be the floppy disk, hard disk, 
magnetic tape, optical disk, optical magnetic disk, DVD, or 
a non- volatile memory card, other than CD-ROM. Also, the 
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program may be downloaded into the computer over a network 
such as the Internet. 

Also, the computer may perform the supplied program to 
implement the functions of the above embodiments, but also 
the program in cooperation with the OS (Operating System) or 
other application software to implement their functions. 
Moreover, all or part of the processing of the supplied program 
may be performed in a function extension board or a function 
extension unit of the computer to implement the functions of 
the above embodiments. These programs are also included 
within the embodiments of the invention. 

The data interpolating device as shown in FIGS. 4 and 
14 performs the data interpolation by sequentially inputting 
a plurality of discrete data, but may be employed as a sampling 
function generating device for generating the sampling 
function as indicated by #4 in FIG. 6 by inputting the single 
digital data. 

Also, the data interpolating device may be applied to 
the digital-analog conversion for converting discrete digital 
data to a continuous analog signal, the expansion process for 
restoring the original data by interpolation from the data 
discretized at the time of compression, and the enlargement 
process of image data, for example. 

In this case, each of the above embodiments only 
illustrates a specific example for practicing the invention, 
but a technical scope of the invention may not be limitatively 
construed by the above embodiments. That is, this invention 
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may be practiced in various ways or forms without departing 
from the scope or spirit of the invention. 

As described above , with this invention , when the discrete 
data are sequentially inputted, interpolation data is 
determined by performing weighted addition of discrete data 
at a sampling point of notice and discrete data at sampling 
points around the sampling point of notice according to the 
value of a digital basic function, whereby the sampling 
function of finite supports dif f erentiable once or more times 
over the whole range is obtained only by very simple 
arithmetical operations, and smooth data interpolation based 
on the sampling function is made by simple arithmetical 
operations. Accordingly, it is possible to not only obtain 
excellent interpolation data without phase distortion or 
truncation error of the low pass filter, but also shorten the 
operation time and simplify the operation circuit. 

[Industrial Applicability] 

This invention allows smooth data interpolation to be 
made based on the sampling function of finite supports 
dif f erentiable once or more times over the whole range by very 
simple arithmetical operations, whereby it is possible to 
shorten the operation time and simplify the operation circuit . 



